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ABSTRACT 

The ROSAT Deep Cluster Survey (RDCS) has provided a new large deep sample of X-ray selected 
galaxy clusters. Observables such as the flux number counts n(S), the redshift distribution n{z) and 
the X-ray luminosity function (XLF) over a large redshift baseline (z^S 0.8) are used here in order to 
constrain cosmological models. Our analysis is based on the Press-Schechter approach, whose reliability 
is tested against N-body simulations. Following a phenomenological approach, no assumption is made 
a priori on the relation between cluster masses and observed A-ray luminosities. As a first step, we 
use the local XLF from RDCS, along with the high-luminosity extension provided by the XLF from the 
Brightest Cluster Survey, in order to constrain the amplitude of the power spectrum, erg, and the shape 
of the local luminosity-temperature, L bo i~T, relation. We obtain a s = (0.58 ±0.06) x f)-°- 47 + - 16 °o f or 
flat models (Q A = 1 - fi ) and a 8 = (0.58 ± 0.06) x ft °- 53 +°- 2 ™° for open models (fi A = 0) at 90% 
confidence level, almost independent of the L^-T shape. The density parameter f2 an d the evolution 
of the L bol -T relation are constrained by the RDCS XLF at z > and the EMSS XLF at z = 0.33, and 
by the RDCS n(S) and n(z) distributions. By modelling the evolution for the amplitude of the Lb a i~T 
relation as (1 + z) A , an flo = 1 model can be accommodated for the evolution of the XLF with 1^5 A^, 3 
at 90% confidence level, while fio = 0.4^° 2 an d 0.6 are implied by a non-evolving L^oi-T (A = 0) 
for open and flat models, respectively. 

Subject headings: Cosmology: theory - dark matter - galaxies: clusters: general - X-rays: galaxies 



1 Introduction 

Galaxy clusters are crucial probes for models describing 
the formation and evolution of cosmic structures. In stan- 
dard scenarios, clusters form at high peaks of the primor- 
dial density field (e.g., Kaiser 1984; Bardeen et al. 1986). 
Therefore, both the statistics of their large-scale distribu- 
tion and their abundance are highly sensitive to the nature 
of the underlying dark matter density field. Furthermore, 
their typical scale, ~ 10/i _1 Mpc (h is the Hubble con- 
stant in units of 100kms _1 Mpc _1 ), relates to fluctuation 
modes which have still to enter, or are just approaching, 
the non-linear stage of gravitational evolution. Then, al- 
though their internal gravitational and gas dynamics are 
rather complex, a statistical description of global cluster 
properties can be provided by resorting to linear theory or 
perturbative approaches. By following the redshift evolu- 
tion of clusters, we have a valuable method to trace the 
global dynamics of the Universe and, therefore, to deter- 
mine its geometry. 

In this context, the cluster abundance at a given mass 
has long been recognized as a stringent test for cosmolog- 
ical models. Typical rich clusters have masses of about 
5 x 10 14 /i~ 1 A/ Q , which is quite similar to the average mass 
within a sphere of ~ 8/1 Mpc radius in the unperturbed 
universe. Therefore, the local abundance of clusters is ex- 
pected to place a constraint on us, the r.m.s. fluctuation 
on the 8 /i -1 Mpc scale. Analytical arguments based on the 



approach devised by Press & Schechter (1974) show that 
the cluster abundance is highly sensitive to a$ for a given 
value of the density parameter £Iq (e.g., Frenk et al. 1990; 
Bahcall & Cen 1992; Lilje 1992; White, Efstathiou & Frenk 
1993; Viana & Liddle 1996). Once a model is tuned so as 
to predict the correct abundance of local (z^ 0.1) clus- 
ters, its evolution will mainly depend on f2o (e.g., Oukbir 
& Blanchard 1992). Hence, the possibility of tracing the 
evolution of the cluster abundance with redshift should 
provide both the value of the density parameter and the 
amplitude of the fluctuations at the cluster scale. 

It is worth stressing, however, that by 'abundance' one 
means the comoving volume density of galaxy clusters per 
unit mass, which cannot be directly observed. Therefore, 
these arguments always assume a relation of some observ- 
able cluster property to the total mass. In this respect, 
X-ray selected samples have been proved to be extremely 
useful, essentially due to the simplicity of their selection 
functions. 

The availability of new observational data for high-z 
clusters has recently stimulated a flurry of activity in this 
direction. For instance, Carlberg et al. (1997) and Fan, 
Bahcall & Cen (1997) have used results on the internal ve- 
locity dispersion a v of X-ray selected clusters in the CNOC 
survey (Carlberg et al. 1996) and concluded that low- 
density models with as — 0.7-0.9 and Oq — 0.3-0.5 are 
preferred (see, however, Gross et al. 1998). 
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Observations of cluster X-ray temperatures offer an in- 
dependent and powerful means to characterize the evolu- 
tion of the cluster abundance since cluster temperatures 
are directly connected to cluster masses. Eke, Cole & 
Frenk (1996) have pointed out that a measurement of the 
cluster temperature distribution at z ~ 0.5 would discrim- 
inate at a high confidence level between a critical den- 
sity model and a low-density model with f^o ~ 0.3. First 
attempts in this direction have been pursued by Henry 
(1998) and Eke et al. (1998) who used available data 
on the X-ray temperature function for z^> 0.4 and found 
us — 0.5-0.8 and flo — 0.3-0.7. Based on a similar data 
set, a different conclusion has been however reached by 
Viana & Liddle (1998), who argued that a critical-density 
Universe is still viable as far as the evolution of the clus- 
ter temperature function is concerned. More robust re- 
sults will require substantially larger compilations of clus- 
ter temperatures, which will be a time consuming process 
even with the next generation of X-ray satellites. 

An alternative way to trace the evolution of the clus- 
ter abundance is to rely on the luminosity distribution 
of X-ray flux-limited cluster samples (e.g., Henry et al. 
1992; Oukbir & Blanchard 1992; Bartlett & Silk 1993; Co- 
lafrancesco & Vittorio 1994; Reichart et al. 1998). The 
advantage of this approach lies in the simplicity of the 
measurement and the availability of large samples, with 
well understood selection biases. A good understanding 
of hydrodynamical processes in the intra-cluster medium 
(ICM), as well as feedback mechanisms from stellar energy 
release, are however needed in order to relate the X-ray 
luminosity to the cluster mass. 

Following the original samples from the Einstein 
Medium Sensitivity Survey (EMSS; Gioia et al. 1990), 
the ROSAT-PSPC has recently given a strong impulse in 
this area providing both large solid angle, high flux limit 
samples (e.g. the Bright Cluster Sample (BCS) by Ebel- 
ing et al. 1997), and more distant cluster samples from 
deeper, small solid angle surveys (the ROSAT Deep Clus- 
ter Survey (RDCS) by Rosati et al. 1995, 1998; WARPS 
by Scharf et al. 1997; SHARC-S by Burke et al. 1997 and 
Collins et al. 1997; the sample based on ROSAT PSPC by 
Vikhlinin et al. 1998). 

In this paper, we use, for the first time, the complete 
amount of information contained in the RDCS, namely 
the evolving X-ray luminosity function (XLF), the deep- 
est flux number counts and the redshift distribution, along 
with a knowledge of the RDCS selection function. In par- 
ticular, in order to normalize the models at z ~ we will 
use the local XLF from the RDCS and the high-luminosity 
extension from BCS. Then, we will use the XLF estimates 
from RDCS at the median redshifts z ~ 0.3 and z ~ 0.6, 
the flux-number counts n(S) and the redshift distribution 
n(z) from the RDCS to constrain the evolution. We note 
that the number counts and redshift distribution represent 
the projection of the ^-dependent XLF along the redshift 
and flux, respectively. The number counts alone from the 
RDCS have been already used by different authors to place 
constraints on cosmological models (see Kitayama & Suto 



1997; Mathiesen & Evrard 1998 for a parametric approach; 
see Cavaliere, Menci & Tozzi 1998 for an approach based 
on a physical models for the ICM). 

Instead of relying on specific physical models to describe 
the ICM and its evolution, we prefer here to adopt a phe- 
nomenological approach. Using a parameteric expression 
for the relation between cluster masses and X-ray lumi- 
nosities, we fit the corresponding parameters, along with 
those describing the cosmological models, to match the 
RDCS data. 

The two principal issues addressed in this paper are to: 

(1) derive robust constraints on the amplitude of fluctua- 
tions on the cluster mass scale using the local XLF and; 

(2) derive robust constraints on Q,q and the L^oi-T rela- 
tion from a flux-limited sample of clusters with z^ 1, like 
RDCS. 

The structure of the paper is as follows. In Section 2 
we give a brief description of the data, while in Section 3 
we present the method of analysis. We review the Press- 
Schechter approach for the mass function and compare its 
predictions with a set of P3M N-body simulations in or- 
der to assess its reliability. Then, we describe in detail 
our procedure to convert cluster masses into luminosities. 
This will be accomplished in three steps: (i) converting 
virial masses into temperatures; (ii) converting temper- 
atures into bolometric luminosities; (Hi) correcting lumi- 
nosities from the bolometric to the soft ROSAT [0.5-2.0] 
keV band. In Section 4 we present the results of this anal- 
ysis from the local XLF to constrain the shape and the 
amplitude a% of COBE-normalized power-spectra as well 
as the local luminosity-temperature (Lboi-T) relation. In 
Section 5 we use the evolving XLF, number counts and 
redshift distribution to constrain the density parameter 
fio and the evolution of Lboi-T. A brief discussion of the 
results and the main conclusions are given in Section 6. 

2 The RDCS sample 

The RDCS compiled a large, X-ray flux limited sample 
of galaxy clusters, selected on the basis of X-ray prop- 
erties alone, via a serendipitous search in ROSAT-PSPC 
deep pointed observations (Rosati et al. 1995). The depth 
and solid angle of the survey were chosen to probe an 
adequate range of X-ray luminosities over a large red- 
shift baseline. Over 160 cluster candidates were selected 
down to the flux limit of S = 1 x 10~ 14 crgs _1 cm -2 , 
over an area of 50 deg 2 , by utilizing a wavelet-based de- 
tection technique. The latter is particularly efficient in 
discriminating between point-like and extended, low sur- 
face brightness sources. The completeness and degree of 
contamination of the catalogue are well understood above 
S = 2 x 10 _14 ergs~ 1 cm -2 , but still uncertain at lower 
fluxes where the optical identification program is not yet 
complete. The sky coverage of the survey, a crucial in- 
gredient in the interpretation of deep cluster surveys, has 
also been extensively studied (Rosati et al. 1998; RDNG 
hereafter). Cluster redshifts have been secured for more 
than 100 clusters/groups to date using NO AO and ESO 
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telescopes. In this paper, we use the flux limited subsam- 
ple of spectroscopically identified clusters used by RDNG 
to derive the cluster XLF. This comprises 70 clusters with 
fluxes S > 4 x ICT^ergs^cm- 2 and 0.05< z< 0.85. The 
whole RDCS sample is used to compare model predictions 
with the observed cluster number counts. An additional 
determination is included in the present analysis at the 
faintest bin (S < 2 x I0~ 14 ergs _1 cm~ 2 ) which extends 
the Log N-Log S presented by RDNG. We note that 15% 
of the cluster candidates, mostly at the faint end, still re- 
main to be optically identified. 

3 The Theoretical framework 

3. 1 The recipe for the cluster mass function 

3.1.1 The power spectrum model 

We write the linear power spectrum of density fluctua- 
tions as P(k) oc k npr T 2 (k). For the transfer function T(fc), 
we assume the Gamma-model 

ln(l + 2.34g) 
T[q) = 2.34 g X 

[l + 3.89g + (16.1g) 2 + (5.46g) 3 + (6.71g) 4 ] ~ 1/4 (1) 

where q = k/hT, being T the shape parameter. For the 
class of CDM models, it is T ~ floh (Bardeen et al. 1986), 
while in general T can be viewed as a free parameter, to 
be fitted to observational constraints. For instance, for the 
fio = 1) case r ~ 0.2 can be obtained either in the frame- 
work of Cold+Hot DM models with a suitable choice of 
the massive neutrino fraction (e.g., Primack 1997), or in 
the framework of rCDM models, where the CDM shape of 
P(k) is modified by the decay of massive neutrinos (White, 
Gclmini & Silk 1995). As for the primordial spectral in- 
dex, n pr , we will mainly concentrate in the following on 
the Harrison-Zeldovich = 1 , although we will also 

comment about the effect of varying n pr around this value. 
As usual, the amplitude of P(k) will be expressed in terms 
of erg and is fixed by the four-year COBE normalization 
recipe, as provided by Bunn & White (1997) and Hu & 
White (1997) for the Q A = 1 - O and the S1 A = cases, 
respectively. This normalization determines a one-to-one 
correspondence between a% and the shape parameter T for 
a fixed value of n pr . 

3.1.2 The Press-Schechter approach 

According to the PS formalism, clusters at a given red- 
shift, z, are identified with those halos that are just viri- 
alizing. The comoving number density of such structures 
in the mass range [M, Al + dM] reads 
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Here p is the present day average matter density and 
5 c {z) is the linear-theory overdensity extrapolated at the 
present time for a uniform spherical fluctuation collaps- 
ing at redshift z. This quantity conveys information 



about the dynamics of fluctuation evolution in a generic 
Friedmann background. It is convenient to express it as 
S c {z) = 5 a {z) [D(0)/D(z)\, where 



D(z) = -n E(z) / — —^dz' 
y ' 2 y ' J z E(z') 3 



(3) 



is the linear fluctuation growth factor. In the above ex- 
pression, E(z) = [{l + z) 3 n + {l + z) 2 n R + n A }^ 2 , where 
= A/3H 2 and flu = 1 — O — (see, e.g., Peebles 
1993). 

The quantity 5q(z) has a weak dependence on f2o for 
both = and SI/? = 0. In the following we will 
adopt for So(z) the expression provided by Kitayama & 
Suto (1996). For a critical-density Universe it is 8 c {z) = 
1.686(1 + 2). 

The rms density fluctuation at the mass scale M, aM, 
is connected to the fluctuation power spectrum according 
to 



' M 
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— / dkk 2 P(k)W 2 (kR) 
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(4) 



Here W(x) is the Fourier representation of the window 
function, which describes the shape of the volume from 
which the collapsing object is accreting matter. The 
comoving fluctuation size R is connected to the mass 
scale M as R = (3M/47rp) 1//3 for the top-hat window, 
W(x) = 3(sinx — x cos x) /x 3 , that we adopt in the follow- 
ing. 

3.1.3 Comparison with N-body simulations 

The reliability of the PS formula has been tested against 
N-body simulations by several authors. As a general re- 
sult, it turns out that eq. (Q) provides an overall satis- 
factory description of the N-body mass function around 
the non-linear mass scale M* (i.e., the mass at which 
S c (z)/aM t — 1), which controls the position of the expo- 
nential cut-off (e.g., Lacey & Cole 1993; Eke ct al. 1996). 

In order to make our own comparison, we ran a set of 
N-body simulations for three different models. Simula- 
tions are run for: (a) an f^o = 1 model with n pr — 0.8 and 
r = 0.35; (b) a flat low-density model with fio = 0.4 and 
T = 0.22; (c) an open model with f2 = 0.6 and V = 0.25. 
Model parameters are described in Table 1. They have 
been chosen in such a way that the corresponding power 
spectra are consistent with the APM galaxy P{k) shape 
(e.g., Baugh & Efstathiou 1993), the local cluster number 
density (e.g., Eke et al. 1996; Girardi et al. 1998) and 
the CMB anisotropics from the 4-year COBE data (e.g., 
Gorski et al. 1996). 

Each simulation is run within a box of C = 250 /i _1 Mpc 
a side using the adaptive particle-particle-particle-mesh 
(AP3M) code, kindly provided by Couchman (1991). The 
code follows the trajectories of 128 3 particles, with a 
comoving equivalent Plummer softening parameter of ~ 
100 h' 1 kpc. Therefore, both the mass resolution (m par t — 
2.06 x 10 12 f2o /i _1 M Q for the mass of each particle) and the 
dynamical resolution are adequate to describe halo masses 
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down to poor clusters (< 10 14 /i -1 M©). The initial red- 
shift Zi at which simulations are started has been fixed so 
that (7 = 0.2 for the r.m.s. fluctuation amplitude on the 
grid at that time. The integration variable is chosen to 
be p = a 3 / 2 , where a = (1 + z)^ 1 is the expansion fac- 
tor. The number of time-steps has been determined by 
fixing Ap = 0.14. This ensured an energy conservation of 
about 3% at the final step of each simulation. For each 
model, we run five different realizations in order to have 
enough statistics to reliably estimate the mass function in 
the high-mass tail. 

TABLE 1 

Parameters of the simulation models 



Model 


n 


/ lp r 


r 


f8 




Time— steps 


EdS 


1.0 


0.8 


0.35 


0.56 


10 


240 


Flat 


0.4 


1.0 


0.22 


0.87 


16 


180 


Open 


0.6 


1.0 


0.25 


0.67 


15 


440 




Fig. 1.— Comparison between the cumulative cluster mass 
function n(> M) produced by N-body simulations (his- 
tograms) and the Press-Schechter formula (smooth curves) at 
different redshifts. Details about the parameters of the simu- 
lated models are provided in Table 1. The horizontal dotted 
line indicates the shot-noise level, i.e. the cluster number den- 
sity associated with one single cluster in the volume of the five 
simulation boxes. The vertical dotted lines indicates the min- 
imum fiducial mass for a cluster to be reliably identified, i.e. 
corresponding to a cluster resolved with 20 particles. 

Clusters in each simulation are identified by following a 
two-step procedure. Firstly, we apply a friend-of-friends 
(FOF) algorithm with linking parameter b = 0.2 in order 
to identify cluster candidates. For each FOF group, we 
compute the center of mass of the particles that belong 
to it. Therefore, starting from each center, draw a sphere 
whose radius encompasses an average density equal to the 



virial density: p V ir(z) — A V i r {z)p(z) where p(z) is the 
cosmic average density at rcdshift z. It is A V i r = 18w 2 for 
Oo = 1, whereas we take the expressions provided by Ki- 
tayama & Suto (1996) for the < 1 cases. The center of 
mass of all the particles within such a sphere is then com- 
puted and used as the new center for the next iteration. 
The procedure is stopped, in general after few iterations, 
once the coordinates of the center of the sphere and its 
mass converge. When two clusters are found to partially 
overlap, the smaller one is discarded from the final list. 

In Figure 1 we plot the cumulative mass function rt(> 
M), defined as the number density of objects with masses 
larger than M. Simulation results (histograms) are com- 
pared to the PS predictions (dashed curves) at four dif- 
ferent redshifts, z = 0,0.21,0.55,1.40. In each panel, the 
horizontal dotted line indicates the shot-noise level - the 
number density corresponding to finding one single cluster 
in all the five boxes. The vertical dotted lines indicate the 
mass of a halo containing 20 particles, smaller halos being 
not properly resolved. 

As a general result, we confirm that the PS approach 
provides a good approximation to the N-body mass func- 
tion in the region around the non-linear mass scale M*, 
quite independently of the model. The ability of the PS 
formula to account for the different degrees of evolution in 
the different cases is remarkable. A similar result about 
the reliability of the PS formula has been also recently 
obtained from the analysis of Hubble-volume simulations 
(White 1998, private communication; see also Colberg et 
al. 1998). Governato et al. (1998) analysed two high- 
resolution simulations, for critical density and low-density 
CDM models, which were aimed at testing the PS formula. 
They found a slightly different redshift dependence of S c , 
which turns into a ~ 6% smaller value at z ~ 0.5 only 
for the f2o = 1 case. Even in this case, such a difference 
only induces variation in as and Qq which are smaller than 
the statistical uncertainties of our analysis (see Section 4, 
below). 

At low masses the PS expression provides a systematic 
overestimate of the halo number density (see also Bryan 
& Norman 1998; Gross et al. 1998). However, given the 
flux-limit of the RDCS sample, such small masses enter 
only at the lowest redshifts, whereas we will mostly use 
the observed distributions at high-redshifts (z<, 0.3) to 
constrain our models (we will further comment in the fol- 
lowing on the mass scales associated to RDCS clusters at 
different redshifts). 

3.2 The mass-temperature relation 

According to the spherical collapse model and under the 
assumption of virial equilibrium and isothermal gas distri- 
bution, the mass-temperature relation can be written as 



k R T = 



1.38 



M 



2/3 



10 15 h- 1 M Q/ 
x [n A mr (z)} 1/3 (l + z)keV, 



(5) 
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where 76% of the gas is assumed to be hydrogen (see, e.g., 
Eke et al. 1996). The (3 parameter is defined as the ratio 
of the specific kinetic energy of the collisionless matter to 
the specific thermal energy of the gas, 

/?=^#, (6) 

being fi — 0.59 the mean molecular weight, m p the proton 
mass and o~ v the one-dimensional cluster internal velocity 
dispersion. Observational data on clusters with reliable 
determinations of both X-ray gas temperature and galaxy 
velocity dispersion indicate (3 ~ 1 for 3 keV, (e.g., 
Lubin & Bahcall 1993; Girardi et al. 1996). The calibra- 
tion of the (3 value using numerical simulations have been 
attempted by several authors (see, e.g., Bryan & Norman 
1998, for a summary of numerical results), by fitting either 
the M-T relation of eq. (||) or the T-a v relation of eq. (j^) . 
Differences between the resulting (3 values could be due to 
a non-perfect virialization of the cluster or to a departure 
from the hydrostatic equilibrium in the ICM. Here we will 
refer to results based on the M-T relation, that we will 
use to convert masses into X-ray luminosities. 

Using a Smoothed Particle Hydrodynamics (SPH) sim- 
ulations of six clusters, Evrard (1991) found (3 = 1.23 by 
adopting a mass-weighted temperature definition. Evrard, 
Metzler & Navarro (1996) analyzed an enlarged sample of 
clusters, also including low-fio cases, and found an av- 
erage value of (3 = 1.05. Their results also show some 
dependence on Oo, low-density models preferring a lower 
(3. Bryan & Norman (1998) analyzed a set of cluster sim- 
ulations based on the piecewise-parabolic-mesh (PPM) 
method. Using a luminosity-weighted definition for T, 
they found 1.27 < (3 < 1.33. Pen (1998) used Moving- 
Mesh Hydrodynamic (MMH) cluster simulations for dif- 
ferent cosmological models. By assuming an emission- 
weighted temperature, he found (3 = 1.12 ± 0.04. Re- 
cent simulations of the Santa Barbara Cluster (Frenk et 
al. 1998), based on a variety of numerical techniques, con- 
verge to indicate that (3 ~ 1.15. This value, which also 
falls within the range of previous results, will be adopted 
in the following as the fiducial one. We will also show 
results for [3 = 1. 

3.3 The luminosity-temperature relation 

The observational determination of the relation between 
bolometric luminosity and temperature, Lb i~T, at low 
redshift is becoming more and more accurate as larger and 
higher quality data sets are constructed (e.g., Mushotzky 
1984; Edge & Stewart 1991; David et al. 1993). If we 
model the Lboi~T relation as 

L b oi = L 6 (l + z) A x W 44 h- 2 erg s" 1 , (7) 

then low redshift data for 3 keV indicates L 6 ~ 3 as a 
rather stable result, and a ~ 2.7-3, depending on the sam- 
ple and the data analysis technique. As for the behavior at 
lower temperatures, Ponman et al. (1996) analysed a set of 



ROSAT observations for 22 Hickson's compact groups and 
found indications for a steepening of the Lboi~T relation 
below 1 keV. White, Jones & Forman (1997) analysed a set 
of 207 EINSTEIN clusters and found a ~ 3. Although the 
formal fitting uncertainties are generally small, the scatter 
of data points around the relation (0) is so large as to raise 
the question of whether it represents a good model for the 
observational Lboi~T relation. In the same paper, White et 
al. also confirm the result by Fabian et al. (1994) concern- 
ing the dependence of the scatter in the Lboi~T relation 
on the cooling-flow (CF) mass deposition rate, a smaller 
scatter being found when excluding CF clusters. Similar 
results about a tight Lboi~T relation have been obtained 
also by Arnaud & Evrard (1998), Markevitch (1998) and 
Allen & Fabian (1998) by either considering only clusters 
with a low CF or by correcting for the CF effect. 

As for the evolution of the Lb i~T relation, Mushotzky 
& Scharf (1997) compared results from a sample of ASCA 
temperatures at z > 0.14 with the low-redshift sample by 
David et al. (1993). They found that data out to z ~ 0.4 
are consistent with no evolution (i.e., A ~ 0), although 
within rather large uncertainties. Henry (1997) deter- 
mined the luminosity-temperature relation for a sample of 
EMSS clusters at a median redshift z = 0.32 using ASCA 
temperature measurements. The decrement he found in 
the amplitude of the Lboi~T relation with respect to that 
by David et al. (1993) implies a marginally positive evo- 
lution with A = 0.36 ± 0.32. Sadat, Blanchard & Oukbir 
(1998) analysed a compilation of z > clusters taken from 
different authors and also found a mildly positive evolving 
L bol -T with 0< A< 1. 

On the theoretical side, the first attempt to describe the 
Lboi~T relation and its evolution has been developed by 
Kaiser (1986) for an Einstein-de Sitter (EdS) model. This 
model, based on the assumptions of self-similar gas evolu- 
tion, can be extended to generic cosmologies and predicts 
Lboi oc T 2 A*^(z). The resulting slope of the local Lboi~ 
T is definitely shallower than the observed one. On the 
other hand, comparisons between the self-similar scaling 
and results from hydrodynamical simulations have shown 
a general good agreement (e.g., Eke, Navarro & Frenk 
1998, and references therein), especially when the effects 
of finite numerical resolution are taken into account (e.g., 
Bryan & Norman 1998). Therefore, the basic discrepancy 
between the predicted Lboi~T relation and the observed 
one, Lboi oc T a with a ~ 3, must be considered as an 
open problem whose solution calls for additional physics, 
like the preheating generated by non gravitational sources 
(e.g., Cavaliere, Menci & Tozzi 1997; see, however, Metzler 
& Evrard 1997). 

As for the z-dependence of the Lb i~T relation, Kaiser 
(1986) predicted an induced evolution in the XLF, for ac- 
ceptable scale-free power spectra, which is opposite to the 
mild negative evolution found later in the data (e.g., Gioia 
et al. 1990). This led Kaiser (1991) and Evrard & Henry 
(1991) to assume that an initial entropy was imprinted in 
the ICM. Afterwards, this model has been generalized by 
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Bower (1997) to include the general redshift evolution for 
the minimal entropy level. As a result, a range of evolu- 
tionary patterns can be derived for Lboi~T and, therefore, 
for the XLF. 

In our analysis we prefer to adopt a conservative ap- 
proach and, instead of assuming a unique shape and evo- 
lution for the Lt> i-T relation, we fit the corresponding 
parameters a and A to the observational data. The am- 
plitude of the Lb i-T relation is taken to be L§ = 2.9 which 
represents the best-fitting value to the data by White et 
al. (1997). 

In order to convert bolometric luminosities into the soft 
ROSAT band [0.5, 2.0] keV we need to introduce the ap- 
propriate correction. We perform it by using a Raymond- 
Smith code and assuming an overall ICM metallicity of 
0.3 times the solar value, as actually observed. The ra- 
tio between the bolometric and the finite-band luminos- 
ity, Lboi/L\o,B— 2.0] ; is plotted in Figure 2. In this plot we 
also show the differences respect to a pure bremsstrahlung 
spectrum with a power-law approximation for the Gaunt 
factor, g(E/k B T) oc (E/k B T)-^. Although this simplified 
model for 7 = 0.3 provides a reasonably accurate bolomet- 
ric correction for rich clusters, it becomes inadequate for 
T^S 2 keV where the effect of metal emission lines in the 
cluster spectra start playing a non-negligible role, increas- 
ing the emissivity above the Bremsstrahlung prediction 
(Sarazin 1988). 




2 4 6 8 10 

T (keV) 

Fig. 2. — The bolometric correction, as computed from 
a Raymond-Smith code with ICM metallicity of 0.3^© 
(solid curve), compared with the prediction of the pure 
brcmssthalung for different power-law approximations for the 
Gaunt factor (short-dashed, dotted and long-dashed curves). 

Therefore, given the above recipe to convert masses into 
luminosities, each model is specified by four parameters, 
namely a^^VL^^a and A. The first two parameters spec- 
ify the cosmological scenario, while the others are linked 



to the thermodynamics of the ICM. We follow a two-step 
procedure to place constraints on such parameters: 

(a) for each value of we constrain T, or equivalently 

erg, and a with the local XLF; 

(b) we use XLF at z > 0.3, number counts and red- 

shift distribution from RDCS to constrain Slo and 
A, which specify the fluctuation growth and the ICM 
evolution, respectively. 

The amplitude of the Lb i-T relation is taken to be 
L$ = 2.9 which represents the best-fitting value to the 
data by White et al. (1997) and it is also consistent with 
the data by Arnaud & Evrard (1998) for clusters with- 
out cooling flows. The remaining parameters to be deter- 
mined are A, which is connected to the thermodynamics 
of the ICM, and f^o (along with which determines the 
growth rate of fluctuations. The results will be presented 
as constraints on the Qq-A parameter space. 
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Fig. 3.— The local XLF from RCDS (filled circles) and from 
BCS (open circles) are compared to model predictions. Left 
and right panels are for an Qo = 1 and a flat fio = 0.3 model re- 
spectively. Upper panels show the effect of changing the shape 
parameter V and, therefore, the as amplitude. Lower panels 
show the effect of changing the shape of the Lboi~T relation. 

4 Results from the local XLF 

The XLF in the soft ROSAT energy band [0.5-2.0] keV 
is related to the PS mass distribution according to 

: frUr MM)dM 

where for simplicity L represents the luminosity in the cor- 
responding band. In their paper, RDNG provided the XLF 
for the RDCS sample at different redshifts for luminosities 
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already converted to the local rest-frame band by apply- 
ing the appropriate X-correction. We refer the reader to 
this paper for more details about the XLF computation. 

In Figure 3 we show an example of how the XLF can be 
used to place constraints on the model parameter space. 
We plot the local XLF from the RDCS (filled circles) and 
the BCS at L x > lO^ergs" 1 (Ebeling et al. 1997; open 
circles) in order to cover the bright end of the XLF, not 
probed by the RDCS. We also show different model pre- 
dictions to emphasize the effect of changing the spectrum 
amplitude as (upper panels) and the shape of the Lm-T 
relation a (lower panels) for a critical density model (left 
panels) and a flat low-density model with SIq = 0.3 (right 
panels). As expected, the effect of changing as in the XLF 
is similar to that in the mass function: a rapid change in 
the exponential cut-off with a smaller and smaller effect 
in the low-luminosity tail. A steepening of the L^oi-T 
relation (i.e., a larger a) corresponds to a larger luminos- 
ity range for a fixed temperature range. As a result, the 
corresponding 4>(L) becomes shallower. 

Constraints on the amplitude of the power spectrum 
(or, equivalently, the shape parameter T) and the slope of 
the local L^oi-T relation are placed by adopting a two- 
parameter ^-minimization procedure. The estimate of 
the x 2 between model predictions and data is performed 
by assuming a log-normal distribution of observational un- 
certainties in (£(£). The probability for a model to be 
accepted is then computed as the probability that the 
data come from the parent model distribution with Gaus- 
sian random variations in logarithmic units given by the 
size of the error-bars. Note that the BCS and the local 
RDCS XLFs have two slightly different median redshifts, 
(z) ~ 0.1 and (z) ~ 0.17, respectively. This leads to a 
~ 6% ambiguity when constraining the amplitude of the 
power spectrum and, therefore, as- 

We show in Figure 4 the results from the local XLF. 
For both the fin = and fl\ = cases (left and right 
panels, respectively) we plot the flo dependence for the 
shape parameter T, the corresponding as from the COBE 
normalization and the slope a of the Lboi—T relation. The 
shaded regions represent the 90% c.l. for model rejection. 
The 90% c.l. for each one of the two fitting parameters 
has been determined by fixing the other parameter at its 
best-fitting value (i.e., minimum x 2 ). The dashed lines 
indicate the minimum-x 2 parameters if (3 = 1 is assumed 
instead of the fiducial value (3 = 1.15. 

As for as and T, we stress that these two quantities 
are connected by a one-to-one relation, once the power- 
spectrum is COB-E-normalized. As we show in Fig. 4, 
they are determined with rather small uncertainties, thus 
confirming that the cluster abundance provides a strin- 
gent constraint on the amplitude of mass fluctuations at 
the cluster scale. The as~flo relation plotted in Fig. 4 can 
be analytically fitted as: 

a s = (0.58 ± 0.02 ± 0.04) x f7- a47+016n ° ; 

£l\ = 1 — CIq 
a s = (0.58 ± 0.02 ± 0.04) x f7 °' 53+0 - 27n ° ; 



n A = o 



(9) 



The first uncertainty corresponds to the 90% c.l. from 
the x 2 minimization, whereas the second error reflects the 
lack of a common median redshift of the BCS and the local 
RDCS XLFs. Systematic uncertainties due to variations 
of S c with respect to the canonical spherical-collapse value 
are not included in eq. ([)]). However, both our analysis 
of iV-body simulations presented and the most recent re- 
sults by other authors (e.g., Governato et al. 1998; White 
S.D.M., private communication), converge to indicate that 
such uncertainties are smaller than the errors quoted in 
eq.(|). 
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Fig. 4. — Constraints on the model parameters from the lo- 
cal XLF, by combining the RDCS and the high-luminosity 
(L[ .5_2] > 10 44 ergs _1 ) BCS data. Left and right panels refer 
to flat (Qr — 0) and open (Qa = 0) cases, respectively. The 
shaded areas indicate the 90% c.l. as deduced from the \ 2 ~ 
minimization procedure, by assuming (3 — 1.15 (see text). The 
dashed curves indicates the best fitting parameters for = 1. 

The above scaling for as translates into an increasing 
trend for T if — and a decreasing trend for Ha = 0. 
Such a behavior is due to the different prescription for 
COBE normalization of open and flat models. For a fixed 
shape of the power spectrum, flat models (e.g., Bunn & 
White 1997) require values of as which are relatively larger 
and larger than those required by open models (e.g., Hu 
& White 1997) as f^o decreases. Therefore, in order to 
compensate for this effects, open models select shallower 
spectra (i.e., larger Ts) at small Slo- We also note that 
taking (3 = 1 decreases the central value in eq.(^) to 0.54. 
This result agrees with that found by Eke et al. (1996) 
who in fact assumed (3 = 1, but used the local X-ray clus- 
ter temperature function. Furthermore, any variation of 
5 C with respect to its canonical value turns into a propor- 
tional change of as ■ 
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Weaker constraints are instead obtained for the shape 
of the Lboi-T relation. For open model we find 3^5 a^S 4, 
roughly independent of f^o, while flat models favor some- 
what smaller values at low ^Iq. Such a difference is re- 
quired in order to compensate for the marginally larger erg, 
for flat models at small fio, corresponding to a shallower 
XLF. The acceptable a values cover the range of current 
observational determinations of the Lboi~T relation. 

4-1 Effect of tilting the spectrum 

The results so far obtained from the local XLF assume a 
scale-invariant primordial spectrum. On the other hand, 
for a fixed Hq, the local XLF depends to a good approx- 
imation only on a^. Therefore, for n pr ^ 1, the shape 
parameter T should be varied so as to leave as unchanged. 
In this respect, all the constraints on A and Qq that we 
will provide in the following are essentially independent of 
the assumption of a Harrison-Zeldovich primordial spec- 
trum. In Figure 5 we show how T must be varied with 
n pr , for different f2o, so as to keep <rg fixed at the best- 
fitting value reported in Fig. 4. For flat models we show 
the — 0.1, 0.3, 0.5, 1 cases, from lower to upper curves, 
while fio = 0.2, 0.3, 0.5, 1 cases are shown for open models, 
from upper to lower curves (a vanishing contribution from 
tensor mode fluctuations to the CMB anisotropy is always 
assumed). For instance, if f2o = L changing n pr = 1 to 
0.8 implies that the shape parameter must be increased 
from r ~ 0.24 to ~ 0.36 in order to provide an equally 
good fit to the local XLF. The shaded region indicates 
the 95% c.l. constraint, T = 0.23 - 0.28(1 - l/n pr ) with 
15% uncertainty, obtained by Liddle et al. (1996) from 
the shape of the APM galaxy power spectrum. For flat 
models we find that n pr ^t 1 would require Q ^, 0.5, while 
a lower density parameter implies n pr ^ 1. For open mod- 
els, a blue n pr > 1 spectrum is required by 0.3^ Slo~ 1 
with a rapidly increasing n pr for lower values. 




0.7 0.8 0.9 1 1.1 1.2 0.7 0.B 0.9 1 1.1 1.2 



Fig. 5. — The constraints from the local XLF on the shape pa- 
rameter r when the primordial spectral index n pr takes values 
different from unity. Each solid curve provides the best-fitting 
Ts for different Qo- For flat models (left panel) results are pro- 
vided for fio = 0.1, 0.3, 0.5, 1 from lower to upper curves, while 
fio = 0.2, 0.3, 0.5, 1 are reported for open models (right panel), 
from upper to lower curves. The shaded area represents the 
95% c.l. constraint from the power spectrum of APM galaxies 
(Liddle et al. 1996). 



5 Tracing the redshift evolution 

In this section we will describe how the evolution of the 
XLF, the number counts and the redshift distribution can 
be used to place constraints in the (Flo, A) plane. 

5. 1 The evolution of the XLF 

In Figure 6 we compare models predictions for two rel- 
evant cases with the XLF data from EMSS by Henry et 
al. (1992) in the redshift bin z = [0.3 — 0.6] (median red- 
shift (z) = 0.33) and from the RDCS in the redshift bins 
z = [0.25 - 0.50] and [0.50 - 0.85] (median redshifts of 
(z) = 0.31 and 0.60, respectively). The XLF for different 
redshift bins have been separated by A log 4>(L) = 2 from 
each other for reasons of clarity. The EMSS XLF has been 
converted to the 0.5-2.0 keV band as described by RDNG. 
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Fig. 6.— The high-redshift XLF from EMSS (open circles) and 
from RCDS (filled circles) are compared to model predictions. 
The EMSS XLF refer to the redshift interval z = [0.3 - 0.6], 
while that from RDCS refers to z = [0.25 - 0.50] and z = 
[0.50 — 0.85] respectively. The left and the right panels are for 
a flo = 1 and a flat Q,o = 0.3 model, respectively. For each 
model and at each redshift, different curves refer to different 
evolutions for the L^-T relation. 

In Figure 6 the shape of the power spectrum is chosen 
so as to fit the local XLF (cf. Fig. 4). As expected, a 
positive evolution of the Lboi~T relation, i.e. an increase 
of the parameter A, produces more luminous high-z ob- 
jects, thus an higher XLF. The higher rate of evolution of 
the mass function in the £!o = 1 case requires a positive 
evolution of Lboi-T in order to match the lack of evolu- 
tion of the XLF over a wide luminosity range. For the 
f2o = 0.3 case instead, the PS mass function evolves slowly 
at z^S f (cf. Figure f ) so that a negligible or mild negative 
Lboi-T evolution is required by the RDCS XLF. A some- 
what more negative evolution seems to be required by the 



8 



steepening of the XLF observed from the EMSS data at 
Lx^j 4 x 10 44 ergs _1 , although the effect is statistically 
marginal. 

We note that the absence of this steepening of the bright 
end of the XLF at high redshifts would imply a more pos- 
itive evolution of the L^oi—T relation. For instance, we 
verified that assuming a non evolving XLF over the whole 
Lx range would require A ~ 2 for an f^o = 1 model in 
Fig. 6, and A ~ 0.5 for the Q = 0.3 flat model. 

5.2 The number counts 

We now consider the differential cluster number counts, 
n(S), i.e. the number of objects per steradiant with flux 
in the range [S, S + dS] . This observable can be computed 
starting from eq.(0) according to 



n{S)dS 



c 



r 2 (z) dM 
dz^n[M(S,z);z] — dS 



(10) 

(cf. Kitayama & Suto 1997) where r(z) is the radial co- 
ordinate appearing in the Friedmann-Robertson-Walker 
metric: 



r(z) = 
r(z) = 



dzE^(z) ; n A = l-tt Q 
2 [fl z + (2 - fio) (1 - VI + ihz)] 



ng(i + *) 



0. 



(11) 



The flux S in the ROSAT band is related to the luminosity 
according to 

s =4w' (12) 

where dL{z) = r{z)(l + z) is the luminosity distance 
at redshift z. RDNG provided n(S) for S > 2 x 
10 -14 ergs _1 cm -2 drawn from a sample which is a fac- 
tor two deeper than the one used to compute the XLF. 
Here we consider a further extension of the Log N-Log S 
to the survey flux limit. The large error bar of the faintest 
data point is due to incomplete optical identification at 
these fluxes. 
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Fig. 7. — The effect of the Lboi~T evolution on the differential 
number counts, n(S), for the same cosmological models shown 
in Figure 6. 

In Figure 7 we show the effect of varying the evolution 
of the Lb i-T relation on n(S) for the same models as in 



Fig. 6. Again, a larger value of A turns into an increase of 
the flux coming from clusters at a given redshift, so as to 
increase the number counts. Quite similar to the results 
from the XLF, the critical-density model favors a weak 
positive evolution while the low-density model requires a 
negative evolution for Lbd-T. 

5.3 The redshift distribution 

The differential redshift distribution, n(z), is defined as 
the number of clusters above a given flux limit, in the ef- 
fective survey area of the RDCS, with redshift in the range 
[z, z + dz]. Its expression is derived similarly to eq.(|lO|) and 
reads 



n(z)dz 



r 2 (z) 
E{z) 



Slir, 



dS fsky(S) 



x n[M(S, z); z] dz . 

dS 



(13) 



In the above expression Su m = 4 x 10~ 14 erg s -1 cm~ 2 
is the limiting flux above which the RDCS redshift dis- 
tribution is complete, while f s ky(S) is the effective flux- 
dependent sky-coverage appropriate for the RDCS sample 
(see Fig. 1 in RDNG). The upper panels of Figure 8 show 
the dependence of n(z) on the evolution of the L^-T rela- 
tion for the same two models previously considered. In this 
case, the high-z tail of the redshift distribution is highly 
sensitive to variations of the parameter A. The rather ex- 
tended tail of the RDCS dn(z)/dz rules out a negative A, 
if fio = 1. On the other hand, a flat model with no ~ 0.3 
requires a mildly negative evolutionary index in order not 
to overproduce high-redshift clusters. 
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Fig. 8. — The effect of the L^-T evolution on the n(z) red- 
shift distribution, for the same cosmological models shown in 
Figure 6 (upper panels). The lower panels show for the same 
cases the mass of the smallest cluster that, at each redshift z, 
has a luminosity larger than Ln m (z) — A-Kd 2 L {z)Sum (where 
" erg s _1 cm~ 



S Hm = ix 10- 14 1 



In the lower panels of Fig. 8 we also plot the mass 
corresponding to the flux limit Su m as a function of the 
redshift for the different choices of the L^-T evolution. 
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This quantity corresponds to the smallest mass involved in 
the Press-Schechter computation of n{z). Two aspects of 
this plot should be emphasized. Firstly, assuming a neg- 
ative A, larger masses at high z are included in order to 
keep the flux constant; vice versa, for A > the luminos- 
ity at a fixed mass increases with z and therefore, much 
smaller masses at high redshift are required for a cluster 
to emit above the RDCS flux limit. Secondly, such a min- 
imum mass already at zft 0.1 keeps values at which the 
Press-Schechter formula is succesfully tested against N- 
body simulations (cf. Fig.l). This means that the small- 
mass regime where the PS approach has been shown to 
overpredict the cluster abundance should not affect the 
results of our analysis. 

5.4 Constraints on the A-Qq plane 

In order to estimate the range of A values allowed for 
each model we apply a one-parameter x 2 minimization. 
For each f^o, the values of eg and a are those correspond- 
ing to the best-fit to the local XLF (cf. the previous sec- 
tion). A global summary of the resulting constraints on 
the £Iq-A plane from the three observational constraints 
is shown in Figure 9. For both flat (left panels) and open 
(right panels) models, the shaded areas correspond to the 
90% c.l. for A. As a general result, these findings confirm 
on a more quantitative ground that Slo — 1 models pre- 
fer a positive Lboi~T evolution, while lower A values are 
required by smaller CIq. 
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Fig. 9. — Constraints at the 90% c.l. on the Qo~A plane from 
the z > 0.3 XLF, number counts n(S) and redshift distribution 
n(z) (from upper to lower panels, for both flat and open models 
(left and right panels respectively). The limiting fluxes of the 
samples (in units of ergs _1 cm -2 ) from which the XLF, n(S) 
and n(z) are derived are also indicated. 

As expected, the constraints provided by XLF, n(S) 
and n(z) are almost coincident, although the XLF tends 



to prefer marginally larger values of A. Small differ- 
ences between such constraints may be expected in prin- 
ciple. Indeed, they come from somewhat different sam- 
ples. The XLF constraint includes the EMSS results 
for £[0.5,2] a t ( z ) — 0.33, as well as the RDCS data 
above Su m = 4 x 10~ 14 ergs _1 cm~ 2 . The n(S), in- 
stead, has been computed for RDCS clusters down to 
Sum = 1 x 10~ 14 ergs _1 cm~ 2 , while n(z) is based on the 
same flux-limit as the XLF. 

A tentative analytic fit of the constraints on the 
A plane from the redshift distribution is provided by the 
expressions 



A - 
A 



(3f7 - 2) ± 1 
(4fi - 3) ± 1 



^i? = 0; 
Q A = 0. 



(14) 



A similar result has recently been obtained, only for the 
open models, by Sadat et al. (1998) using the redshift dis- 
tribution of the EMSS sample. As a general result, it is 
evident how an independent determination of the evolu- 
tion of the Lb i-T relation and a deeper understanding of 
the nature of its scatter can turn into a constraint on the 
density parameter. For instance, if a redshift-independent 
Lboi~T relation is assumed (i.e. A = 0), then the evolution 
of the XLF implies n = QA±° Q 3 2 if n A = and fl < 0.6 if 
= 0. These findings are consistent with those obtained 
by Eke et al. (1998) from the analysis of the A -ray tem- 
perature function at zjS 0.4, also based on a non-evolving 
Lboi~T relation. 

As for the physical meaning of the evolution param- 
eter A, detailed models involving both the gravitational 
shock heating and non-gravitational heating (see Cava- 
liere, Menci & Tozzi 1997, Tozzi & Norman in prep.), pre- 
dict a range of A values between and 1.5 at the clus- 
ter mass scale, and slightly lower values, —0.5^5 Afz 1, for 
groups. This result holds if the non-gravitational heat- 
ing rate has a smooth dependence on the redshift, as sug- 
gested by the average star formation rate derived from 
deep galaxy surveys (e.g., Madau et al. 1996). However, 
the duration and the epoch of the heating of the ICM, 
which is not well known at present, critically affects the 
effective value of A, so that a self consistent modeling of 
the feedback processes from galaxies is needed in order 
to constrain A within a narrow range. This justifies our 
present choice to leave A as a parameter free to be fit by 
the evolution of the XLF. 

6 Conclusions 

In this paper we have used observational results from 
the Rosat Deep Cluster Survey by Rosati et al. (1998) to 
place constraints on cosmological models. The unprece- 
dented extension of the RDCS both in redshift (z^ 0.8) 
and fluxes (S> 1 x 10~ 14 ergs^cuT 2 in the [0.5-2.0] keV 
band) makes it, in principle, a suitable baseline over which 
the evolutionary pattern of the cluster abundance and, 
therefore, the density parameter f^o, can be investigated. 
A necessary ingredient for this analysis is the physics of 
the intra-cluster gas, which drives the relation between the 
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distribution of cluster halo masses and the observed distri- 
bution of X-ray luminosities. To this purpose, we adopt 
a parametrical approach in which the shape of the local 
Lboi~T relation, its evolution and the parameters specify- 
ing the cosmological models are all fitted against observa- 
tional data. 

Firstly, we have used the local XLF from RDCS and the 
extension at high luminosities from the XLF of the Bright- 
est Cluster Sample (Ebeling et al. 1997) to fix the r.m.s. 
fluctuation amplitude at the cluster mass scale and the 
shape of the local L^oi-T relation. Secondly, we have used 
the evolving XLF, the flux number counts and the red- 
shift distribution with the aim of constraining the density 
parameter f2o an d the evolution of the L^-T relation. 

The main results of our analysis can be summarized as 
follows. 

(a) A careful comparison between the PS predictions and 

N-body simulations confirms that this analytical ap- 
proach for the distribution of virial halo masses is 
adequate over a rather large mass range around the 
non-linear mass scale. The low mass tail, where 
the PS formula overpredicts the halo abundance, has 
been shown to lie outside the mass range probed by 
RDCS clusters. 

(b) The local XLF data constrain the amplitude of the 

power spectrum according to 

a s = (0.58 ± 0.06) x ^ a47+016n ° ; 

f2j\ = 1 — f2o 
a s = (0.58 ± 0.06) x fi-0-53+o.27fi . 

ft A = 0. (15) 

This result agrees with analyses of the X-ray tem- 
perature function (e.g., Viana & Liddle 1996; Eke et 
al. 1996; Markevitch 1998), the optical virial mass 
function (e.g., Girardi et al. 1998) and the X-ray 
cluster number counts alone. 

As for the shape of the local L^oi-T relation, we 
find 3^5 ais 4, quite independent of the cosmological 
model and in general agreement with observational 
results (e.g., David et al. 1993; White et al. 1997). 



(c) From the observed evolution of the XLF we constrain 
the density parameter f2o and the evolution param- 
eter A of the Lb i-T relation. We find that f^o = 
1 models require a positive evolution of the T->bol~ 
T relation with 1< A< 3. A non-evolving Lh r—T 
(A = 0), consistent with data at z^S 0.4 (Mushotzky 
& Scharf 1997), implies n = 0.4±g % for open models 
and £!o~ 0.6 for flat models. 

This paper set out to address the two principal issues 
given in the Introduction. In summary, the above results 
lead us to conclude that: (1) Available data on the local 
cluster XLF place robust constraints on as as a function 
of flo [cf. eqs.(^|)], without any a priori assumption about 
the local L^-T relation and; (2) Present uncertainties 
in the evolution of the Lf, ;-T relation do not as yet give 
strong constraints on Qq, even with a high-redshift sample 
as deep as RDCS. 

Samples like RDCS will represent a fundamental basis 
for the study of distant clusters in the years to come. The 
next generation of X-ray satellites and the already avail- 
able large optical telescopes should open the possibility of 
determining cluster masses via X-ray temperature mea- 
surements, virial analysis and gravitational lensing stud- 
ies. Carrying on such observations for an even limited 
number of clusters, extracted from a well defined statisti- 
cal sample, will determine both the evolution of the cluster 
internal dynamics and the value of the cosmological den- 
sity parameter. 
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